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Abstract 

Recently, it has been shown that the ground-state energy of a quantum 
many-body system can be written in terms of cumulants. In this paper 
we show that the energies of excited states can be expressed similarly. 
These representations are suitable for various approximation schemes, 
e.g., projection techniques. The explicit use of cumulants ensures size 
consistency in all approximations. 

The theory is then applied to the computation of the energy bands 
of a semiconductor with diamond structure. As we focus especially on 
the effects of electron correlations, we consider a simplified model and 
re-derive results obtained previously by the variational Local Ansatz 
method. 
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I. INTRODUCTION 



Energy is a size-extensive quantity, i.e., the energy of two well-separated, but 
otherwise identical systems equals twice the value of a single system. Approxima- 
tions which are applied when energies of a system are calculated must preserve this 
property. In the standard diagrammatic approach size consistency is ensured in 
any approximation by considering linked diagrams only |^. However, the diagram- 
matic description makes essential use of Wick's theorem, which is applicable only 
if the dominant part of the Hamiltonian is a one-particle operator. Therefore this 
approach is restricted to weakly correlated systems in which the electron-electron 
interaction is weak and may be treated perturbatively. 

It is known from classical statistical mechanics that size consistency is attained 
by expressing size-extensive quantities in terms of cumulants 0. Kubo defined 
the cumulant of a quantum mechanical expectation value 0. It is denoted by a 
superscript ^ and, with respect to the vector it is given by 

d d 



[1.1] 

=A„=0 



One immediately finds, e.g.. 



1.2) 



(01^1^210)^ = (0|AiA2|0)-(0|A|0)(0|A2|0) . 

Based on the Horn-Weinstein theorem |Q, an infinite series termed "connected 
moment expansion" (CMX) has been derived for the ground-state energy ||5|-0. 
Involving cumulants only (rather than ordinary expectation values), any truncation 
of the series leads to size-consistent estimates for the ground-state energy. 

Here we would like to proceed as in Ref. P, where a different representation of 
the ground-state energy in terms of cumulants has been derived and which we refer 
to for details. The interacting many-body system is described by a Hamiltonian H. 
We assume that the ground state |$o) of ^ part Hq of the Hamiltonian H is known 



Ho\%) = eo\%). (1.3) 

Provided that |$o) has a finite overlap with the ground state of H, the ground-state 
energy Eq of H is given by 

^0 = {H\Q) . (1.4) 

Here we introduced the abbreviation 

{A\B) = (<l>o|At5|$o)'= . (1.5) 

It serves as a form on the vector space (the Liouville space) spanned by the operators 
which are defined on the usual Hilbert space, is a superoperator which acts on 
the Liouville space. It is defined by 

f]=l + lim TTT^i • (1-6) 

Lq is the Liouville operator with respect to Hq. It acts on any operator A according 
to LqA = [Hq,A]^. Q characterizes the transformation of the unperturbed ground 
state |$o) into the ground state of H, |^o). Thus it bears similarities to the wave 
operator of quantum mechanics and we will refer to it as "wave operator" in the 
following. Note, however, that Q appears only in an expectation value, Eq. (|1.5| ), 
and does not yield the wave function |\l/o). 

The explicit use of cumulants in Eq. ( |1.4|) ensures size consistency in any subse- 
quent approximation and can be considered as replacing the afore-mentioned "linked 
diagrams". However, it is important to notice that the cumulant representation of 



the ground-state energy, Eq. ( |1.4| ), does not have the limitations which the dia- 
grammatic approaches often have since it holds for an arbitrary decomposition of 
H into Hq and Hi. In particular, Hq need not be an one-particle Hamiltonian (e.g., 
Hartree-Fock) as in the standard diagrammatic approach. Hence, the present de- 
scription is applicable to strongly correlated systems as well, e.g., to spin systems 
described by the Heisenberg Hamiltonian. 



Another advantage of using Eq. (|1.4| ) is that quite different types of approxi- 
mations can be apphed than in diagrammatic approaches. One example is a parti- 
tioning of the Liouville space, which was pioneered by Lowdin and by the Uppsala 
and Gainsville Quantum Chemistry Groups P-pJ]j. Consider a set of operators Ai,, 
which span a subspace of the Liouville space. Taking into account this "relevant" 
subspace only and, thereby, neglecting the remaining part the following ansatz for 
the wave operator Q ( |1.6|) is made 

Q=l + J2vuA,. (1.7) 

V 

The coefficients are determined by 

{A^\R^) = . (1.8) 

This equation holds generally for all operators A and is easily proven by inserting 
the definition ( |1.6| ) of Q [0. Making the linear ansatz ( |1.7| ) for Q is equivalent to 
applying projection technique to Eq. ( |1.4| ) |^. In particular, if we choose the powers 
of the Hamiltonian as relevant operator set, i.e., A^, = H'^ {u = 1 . . . N) we recover 
the improved CMX as in Refs. [§,0. For a more general ansatz than Eq. ( |1.7| ) and 
other approximation schemes, see Ref. | |12| |. 



For these reasons it would be desirable to express not only the ground-state 
energy, but also the excitation energies in terms of cumulants. The way this has to 
be done can be conjectured by considering the Rayleigh-Schrodinger perturbation 
theory. By expanding Q in terms of Hi and taking the limit z ^ one arrives at 
the perturbation expansion of the ground-state energy: 



Eo = eo + Y.[H- 



n=0 



-^Hi] ) . (1.9) 







The generalization to excited states, which suggests itself, is to take the expectation 
values with respect to an excited state of Hq instead of |$o) 

Hol^i) = e^l^i) (1.10) 



so that the expansion of the energy of an excited state of H becomes 

1 



Ei = e, + Y.[H_ 



?i=0 



'l.ll] 



The subscript i indicates that the expectation values are now taken with respect to 
the excited state 



;i.i2) 



[cf. Eq. (|1.5|) ]. Summing the geometric series in Eq. ( [1 . 1 1| ) , we arrive at the conjec- 
ture 

1 " 



E,, 



H 



1 + hm — , ^ , 
z-{Lq + Hi] 



Hi 



;i.i3) 



where VL is defined as in Eq. (O). Note the interchange of the hmits z ^ {] and 
n — s> oo. 

In the next section, we show that the conjecture ( |1.13| ) is indeed true under cer- 
tain circumstances. In Sec. |TTl| we present a specific application of the formalism by 
calculating the quasiparticle energy bands of a semiconductor in the Bond-Orbital 



Approximation (BOA). In Sec. |^ we compare the results with previous work. Fi- 
nally, we give a brief summary and conclude in Sec. 0. 



II. EXPRESSING EXCITATION ENERGIES IN TERMS OF 

CUMULANTS 



In this section we prove that Eq. ( [L.13|) holds indeed provided |($j|\E'j)p > 1/2. 
Here, |\E'j) denotes an eigenstate of H with eigenvalue E^. Further, let = 
I ($.j|\E'^) p. We consider the following function 



/.(t) = log(<l>.|e-^*^|<l>,) = log^/3«e 



(2.1) 



for t > 0. Note that the derivative of fi{—iX) with respect to A is the function 
introduced by Horn and Weinstein in Ref. [0], which tends to the ground-state 



energy at the limit A ^ oo (provided that ^ 0). However, the function defined 
here need not tend to the ground-state energy a.s t —>■ oo since the arguments of the 
exponentials are imaginary. This is shown next. 

The argument of the logarithm is a complex number and, therefore, we have to 
specify the branch of the logarithm. This is done by requiring that fi{t) is analytic 
in a neighborhood of the positive real t-axis. As we assume that 13? > 1/2, it 
is exp(— 7^ for all t > and this requirement can be fulfilled. We 
continue straightforwardly 

m = log^/^^^e-^*^'^ 



log 



{i)„-itE, 



pre 



3^^ 



PI 



-it{Ef,-Ei) 



-itE, + log I3f + log I 1 + ^ ^ -m,-E,) I ^ 27rm(t) 



P^ 



(2.2) 



P 



The integer function n{t) guarantees that we remain on the correct branch of the 
logarithm. n{t) is fixed by the previous requirement, namely that fi{t) is analytic 
for t > 0. The first two terms of the last line in Eq. ( |2.2|) are analytic. The third 
term is analytic, too, because 



E 



n^^-it{E,,-Ei) 



P 



PI 



1 



1 < 1 



(2.3) 



(remember that /jf"* > 1/2). Therefore n{t) is a constant, which vanishes because 
n{0) = 0. As the argument of the logarithm in Eq. ( ^.2[ ) is bounded the asymptotic 
behavior of fi(t) is 

Mt) 



it 



oo 



(2.4) 



i.e., the energy of the excited state is given by the behavior of fi{t) for large t. 

In order to calculate Ei with the help of Eq. (|2.4|), we consider the half-sided 
Fourier transform of 



dt e^'^m 
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(2.5) 



Since fi(t)/t approaches a constant for t ^ co, the integral is well-defined for 
Im z > 0. From the asymptotic behavior ( p.4| ) we readily obtain the energy of 
the excited state in terms of fi{z): 



Ei = lim z'^fi{z) . 



We finally introduce cumulants by making use of Eq. (|1 . 1| ) 



(2.6) 



/,(t) = log($,|e-*^*|$,) = (<l>,|e 



-iHt 



(2.7) 



This enables us to take the Fourier transform (|2.5|) explicitly and after some algebra 
(cf. Ref. 0) we obtain 

1 



H, 1 



z-{Lo + H,] 
Hence the energy of the exited state is given by 



E,= hm z^f,{z) = {H\n), 

2— >iO+ 



(2i 



(2.9) 



where we again used the definition ( |1.6|) of Q and the abbreviation (|1.12|) . 



III. APPLICATION TO SEMICONDUCTORS 



The formalism presented in the previous section will now be used to calculate 
the electronic excitation energies of a semiconductor. These quasiparticle (quasi- 
hole) energies are defined as the energy difference between an excited state with 
momentum k of the N + 1- {N — 1-) particle system and the ground state of the 
A^-particle system, where N denotes the number of electrons in the charge neutral 
system. 

To be specific, we treat a semiconductor of diamond-type structure. It is de- 
scribed by orthogonalized atomic-like sp^ hybrids and the resulting functions are 
denoted by gi{r_)- The corresponding electron creation (annihilation) operators aj^. 
(a,v,) fulfill the usual anticommutation relations. 



In this basis the Hamiltonian reads as follows: 

H = Hq + Hint 



(3.1) 



i,j,k,l,a,a' 



Here 



= J d\ gl{r) (-^A + ^extdl)) g,{r) (3.2) 

is the bare hopping matrix element between the hybrid orbitals i and j, and Kxt(£) 
is the electrostatic potential set up by the nuclei and the core electrons. 

V,,ui = J d'rd'r' g*(r)g^{r)j-^^gl{r!)gi{r!) (3.3) 

are the interaction matrix elements. 

In the Hartree-Fock (HF) approximation, the effective one-particle Hamiltonian 
is given by 

HuF = ftAc7(^ja (3-4) 
where fij denotes the Fock matrix: 

h = + E [^roU - ('^''1 «L'«/.' I • (3-5) 

1$^^) is the HF ground state which has to be determined self-consistently from 
Eqs. ( |3.4| ) and ( |3.5| ). The HF ground-state energy is given by 



E^' = ($^^ \H, + i/int|$"') + i?NN (3.6) 

where E'nn denotes the interaction energy of the nuclei and the core electrons. 

To proceed further, we introduce bonding and antibonding functions for each 
bond / with corresponding creation operators 
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^ , (3.7) 



The operators a\i^ and aljo- refer to the sp^ hybrids II and 12 which form bond 
I. In the bond-orbital approximation (BOA) [jl3[ the iV-particle HF ground state is 
approximated by occupying all bonding states 

l*^^) = n^l|0). (3.8) 

I,a 

In a more refined calculation the bonding orbitals are constructed from the HF 
orbitals by applying the localization method of Foster and Boys [0] or with the 
help of a localization potential, which acts within the space of the occupied HF 
orbitals only. 

When an electron with momentum k is added to the system the + 1-particle 
HF wave function is given within the BOA by 

l^fL) = ALJ^""^) (3.9) 

where AI^^ creates an electron with momentum k in the conduction band m. The 
operator AI^^^ is decomposed into 

ALa = T.»kmAl)4. (3.10) 

and the akmcr{I) are found from the solution of the HF eigenvalue problem. 

The wave function (|3.9|) yields an energy eigenvalue E^^^. From this we find the 
energy of a quasiparticle in the conduction band 

^HF _ rpHF _ rpHF fi TT\ 



which, according to Koopmans' theorem |]T5|], coincides with an eigenvalue of the 



one-electron Hamiltonian Hup. Similarly, the valence bands are calculated. 

We now want to determine the correlation contributions to these quasiparticle 
energies. They reflect the influence of the residual Hamiltonian 



-f^res — H — i^HF 

We have to include the effect of correlations on both, the A^-particle ground-state 
energy, £'0, and the (excited) + 1-particle energy Ekma- 



The first problem has already been treated in Ref. [|16[. The authors made a 
variational ansatz for the ground-state wave function in order to include correlation 
effects. We want to re-derive their results by directly making use of the cumulant 
representation of the ground-state energy, Eq. (|L^ 



= {H\n) . (3.13) 

It is impossible to determine the wave operator Q exactly. Instead, we make an 
ansatz for Q in the Liouville space taking into account only the most important 
processes which are induced by -ffres- Within the Local Ansatz method they are 
described by the following operator |TB 



S1j = Ij:a'j.'Bj.'AIb,^. (3.14) 

S]j excites electrons in bonds I and J from a bonding into an antibonding state, 
thus changing the charge distribution in these bonds. This reflects the van-der- 
Waals interaction between the bonds / and J, as well as the interaction within one 
bond (/ = J). Our ansatz for the wave operator in Eq. ( p.l3|) reads [see Eqs. ( p..7|) 
and (IL^)]: 

IJ 

The coefficients rfj'] are determined from 

Q = {S]j\Hn) . (3.16) 
Thus we are lead to the same result for the ground-state energy as in Ref. |TB|, which 



we refer to for the further evaluation of the expectation values. 
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This procedure to calculate the ground-state energy is advantageous because it 
applies for excited states, too, as is shown below. 

The energies Ekma which we consider next are, in fact, excited energies of the 
+ 1-particle system. According to Sec. they are given by 

Ekma = {H\fl)tma , (3-17) 

where the expectation values are taken with respect to an excited 1-particle state 
l^^ma)- Again, we make an ansatz for the wave operator Q. In this case, however, 
there are two different correlation effects which we have to take into account. They 
have already been discussed in detail in Ref. [|TB|. The first process is the induced 
dipole-dipole interaction as it was present in the A^-particle ground state above. It 
is described by the operators S^j, Eq. (|3.14| ). The second process is the formation 



of a polarization cloud around the added particle. It is described by 

Sh = lT.^\.'Bj.'AlAj^. (3.18) 

a, a' 

This is seen from 

ShKL) = ^^EAl>B,^^Al\<^n (3.19) 

^ a' 

where the extra electron in bond / induces a transition of an electron in bond J 
from a bonding to an antibonding state resulting in a polarization of bond J. 
For the wave operator Q in Eq. ( |3.171) we make the following ansatz 



= 1 + E Vij{kma)Sh + ^ 7rjj{kmcr)S]j . (3.20) 
ij IJ 

Note that the coefficients describing the dipole-dipole interaction, rju, do not neces- 
sarily equal 77} j. Furthermore, the coefficients T]jj{kma) and 7^ij{kma) may depend 
on the quantum numbers k, m, and a. They are determined from 

= {SUmkma (3.21) 
= (SUmkm. ■ (3.22) 
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This linear system of equations for rijj{kma) and 7rij{kma) leads via Eqs. (|3.20|) and 



(p.lTp to improved energies of the + 1-particle system. Subtracting the ground- 



state energy Eq of the A^-particle system, Eq. ( |3.13| ) with Eq. ( |3.15| ), we arrive at 



quasiparticle energies ekma, which include the main effect of electron correlations. 

Obviously, the case of one added hole can be treated in close analogy. This will 
yield the valence bands. 

IV. DISCUSSION 

In this section we want to compare the results with previous work done in 



Ref. |jT6|. There the same model is treated as here. Starting from a HP calcu- 
lation, the authors make a variational ansatz for both, ground and excited states. 
The correlation operators are the same as in Eqs. (|3.14|) and ( |3.19|) . The authors 
consider the quasiparticle energy directly and, by making use of a linked-cluster 
theorem, they introduce cumulants ensuring size-consistency. Finally, they derive 
equations for the variational parameters rjjj and tt/j by minimizing the energy. The 
final equations look similar to the Eqs. (|3.21|) and (|3.22|) derived here, the main dif- 



ference being that AI^^ is subjected to cumulant ordering there, while it is included 
in the wave function here and, thus, not subjected to cumulant ordering. This dif- 
ference is not surprising: When the quasiparticle energy is calculated directly, only 
processes which are connected to the extra particle contribute, e.g., the blocking of 
ground-state correlations in the neighborhood of the extra particle. This is ensured 
by subjecting A^^^ to cumulant ordering. In this work, however, we calculate the 
energy of the A^- and N ± 1-particle system separately. Therefore, we must also 
include processes which are not connected to the extra particle, e.g., ground-state 
correlations far apart from the extra particle. 

By exploiting the properties of cumulants we may rewrite the expectation values 
in Eqs. ( |3.21| ) and ( |3.22| ) such that they are taken with respect to |$o) only, thus 
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subjecting AI^^ to cumulant ordering. Taking the energy difference between E, 



kmcr 



and Eq, we arrive at essentially the same quasiparticle energies as in Ref. |^ as will 
be shown next. 

We begin with the expectation values which involve S]j: 

= iS]jAl^jHALj + {Sl\H). 

Similarly, 

iSlj\HS^^)kma = ('S'/J^Imo- ~ mcr] 'S'^L^Imcr ) + {^IjWS^i) (4.2) 



since {Al^^\HAl 



kmcr) 



"^mcr- The operators S'Jj are connected to A\^^. Hence, 



HAl 



km<T 



-n '-'KL^kma 



(4.3) 



{Sh\H) kma 

When we insert Eqs. (|]lD-(|]3) into Eq. (|]2T]) we find 

+ Y.^KL{kma)[{S]jAi^^\[H - eZ.]Sl,AL^) + {Sh\HSl,)} u^) 

KL 

+ T.^KLikma){S]jAl^^\HSl,Al^^) . 

KL 

The individual expectation values differ in their behavior as N oo. While the 
cumulants involving A^^^^ turn out to be of order A^~^, the others are 0{N^). To 
leading order in 



= iS]j\H)+j:VKLikma)iS]j\HSl, 



(4.5) 



KL 



Comparing with Eqs. ( p. 151) and ( p.l6| ) we find that riKL{kP^<^) is the same as for 
the A^-particle system (in leading order in A^): 



(4.6) 
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(4.7) 



This result was obtained in Ref. |T^, too. 

Next we consider Eq. (|3.22| ), which determines the coefficients itjj. They describe 
the formation of the polarization cloud around the added particle. Inserting again 
Eqs. we find 

KL 
KL 

The only difference to Ref. is that instead of the exact quasiparticle energy tk-ma 
the Hartree-Fock energy e^^^ enters Eq. {^7j\ ). This is due to a different introduction 
of cumulants and different approximations thereafter: As already mentioned above, 
the energies of the A^- and the ± 1-particle system are calculated separately here, 
whereas the quasiparticle energy is calculated directly in Ref. [|^. 
Finally, the quasiparticle energies are given by 



kmcr 



Eo = iH\n),,^^-iH\n) 



H 



1 + E Vijikma)S]j + E 'r^ij{kma)SJj j 

I J IJ ' kma 



H 



Applying Eq. ( [4. 61 ) this can be rewritten in the form 

euma = eZa + E vf] (4™. |^^7j^L.) + E ^ijitma) (^L. \hSJjAI^^) (4.9) 
IJ IJ 

The same expression for the quasiparticle energy was derived in Ref. [Jl6[. By using 
77! J from a ground-state calculation and determining nij{kma) from Eq. (|4.7|) the 
energy bands for diamond are obtained. The energy gap is found to be only 50% 
of the HF energy gap and the widths of the valence and conduction bands are 



considerably reduced from their HF values |16 . 



V. SUMMARY AND CONCLUSIONS 



In this paper, we showed how excitation energies are expressed in terms of cumu- 
lants. This work generalizes the cumulant representation of the ground-state energy 
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which has been derived recently |^ . Introducing cumulants exphcitly imphes that all 
subsequent approximations are a priori size consistent. We would like to emphasize 
that these representations of ground-state and excitation energies hold irrespectively 
of the particular splitting of the Hamiltonian H into a solvable part Hq and the re- 
maining part {H — Hq)^ which has to be treated approximately. This contrasts the 
standard diagrammatic approaches, in which Hq has to be a one-particle operator. 
Therefore, forming the cumulant can be viewed as considering "linked diagrams" 
only in the case of weakly correlated systems and as their generalization in the case 
of strongly correlated systems for which one-particle wave functions would be a poor 
starting point. 

To give an application of the theory, we calculated the quasiparticle energy bands 
of a semiconductor. We made two separate calculations, one for the ground-state 
energy of the charge-neutral system and one for the (excited) system where one 
electron (hole) with momentum k is added. With the help of the new representa- 
tion of excitation energies it was possible to perform both calculations within the 
same formalism, i.e., to apply the same approximations in both cases. Taking the 
difference of the two energies, we obtained the quasiparticle energy. Finally we dis- 
cussed the connection of the present calculations with a variational treatment of this 



problem [|T^] and essentially reproduced the results obtained there. 
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